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ABSTRACT 



Accurate model stellar fluxes are key for the analysis of observations of individ- 
ual stars or stellar populations. Model spectra differ from real stellar spectra due to 
limitations of the input physical data and adopted simplifications, but can be empir- 
ically calibrated to maximise their resemblance to actual stellar spectra. I describe a 
least-squares procedure of general use and test it on the MILES library. 
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1 INTRODUCTION 

Population synthesis models build on both the theory of stel- 
lar structure and evolution, and stellar atmospheres. The 
spectrum of a given stellar population, characterised by a 
single age and particular chemical composition, is the com- 
posite light from stars with a wide range of mass in pro- 
portions consistent with an adopted initial mass function. 
In the calculations, the spectra of individual stars are sam- 
pled from a library, which is either obtained from radiative 
transfer calculations in model atmospheres (see, e.g., Schi- 
avon, Barbuy & Bruzual 2000; Gonzalez-Delgado et al. 2005; 
Coelho et al. 2007) or from observations (see, e.g., Bruzual 
& Chariot 2003; Le Borgne et al. 2004; Vazdekis et al. 2010). 
Such libraries are also used to classify or parameterise indi- 
vidual stellar spectra. 

The use of libraries of observed spectra has the clear ad- 
vantage that the fundamental building blocks for the models 
correspond to real stars, and bypass all approximations in- 
volved in computing model atmospheres and synthetic spec- 
tra, including the need for accurate atomic and molecular 
data. On the other hand, using real stars has the limita- 
tion that they are sampled from the solar neighbourhood, 
and therefore are limited to what nature provides in this 
small region of the universe, and their distribution in the 
parameter space is far from uniform, sampling some regions 
of interest quite sparsely. 

An interesting option is to combine the best of both 
worlds, using a complete and fairly well-sampled grid of 
model spectra, and tying it to a library of observed spec- 
tra. Such a calibration can be formulated as a least-squares 
problem, where we allow smooth corrections to the model 



grid in order to match as close as possible the available ob- 
servations. In this paper we develop a procedure to this end, 
and apply it to the MILES spectral library described by 
Sanchez-Blazquez et al (2006) and Cenarro et al. (2007). 



2 PROCEDURE 

We consider modelling the differences between synthetic 
spectra of stars and the available observations with a smooth 
function. In particular, we propose to evaluate at each fre- 
quency the ratio between the observed spectra (O a ) and the 
models (M x ), and approximate its dependence on the stel- 
lar parameters by polynomials. We consider spectral energy 
distributions, i.e. spectra that preserve the shape of the con- 
tinuum, although the problem could be formulated as well 
with continuum-normalised spectra. 

We assume that there are only three atmospheric pa- 
rameters characterising the model spectra: the surface tem- 
perature (Toff) and gravity (logg), and the overall metallic- 
ity ([Fe/H0). These parameters are transformed into nor- 
malised quantities for convenience. For T e g we define 



X = (T ff - min(T c s))/(max(T e H) — min(T e s)), 



(X) 



and similar transformations are applied to logy and [Fe/H] 
to define the variables Y and Z, respectively. 

Then our problem is to minimise, for each frequency, 
the merit function 



E 



Os/M s - YjYj a ijk XtYiZ h s 



(2) 
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1 [Fe/H] = \og(N Fe /N H ) - \og(N Fe /N H ) C:J 

2 In the following we will drop the superindex A for simplicity 
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where the index s runs through all the observed stars, and 
the indices i, j, and k, run from to n, to m, and to 
h, respectively. Thus, n, m, and h define the orders of the 
polynomial in X, Y and Z, respectively. The model fluxes 
M s may be calculated directly for each star, but in our case 
they will be obtained by interpolation in a library (see 

To find the minimum, we calculate the derivatives rela- 
tive to each of the parameters of the polynomial (a a ,a 7 ) and 
equate them to zero, arriving at 

E E E °«* E xi +a Yt z^ = ]r ^-x«Yfz:.(3) 

i j k s s 

The indices a, /3, and 7, also run from to n, to m, and to 
h, respectively, leading to a system of (n+l) x (m+1) x (h+1) 
equations. 

Eq.[3]can be written in a more compact way if we map 
the sets of indices ijk and a/?7 into a pair of indices, becom- 
ing 

(n + l)(m+l)(Ji + l)-l 

£ ^V<V = tv, (4) 

j'=o 

1 = 0, f , (n + l)(m + l)(h + 1) - 1, 
where 

y lY = ^x' +a yf^ s fc+ \ (5) 

3 

Cj> = a ijk , (6) 

b *> = Y,it x ° Yfz ' J > (7) 

s 

with 

L/7((m + l)(A + l))J (8) 
3 = j'-[i(m+l)(h+l))/(h+l)\ 

k= j' -i(m+ l)(h+ 1) -j(h+ 1), 
and similarly 

L*7(("» + i)(ft + i))J (9) 

/3= i'- [a(m + l)(A+l))/(ft+l)J 
7= i'-a(m + l)(/i + l)-/3(ft+l), 

where [xo J indicates the integer part of w. 

Eq. 4 can be solved with any of the standard numerical 
methods for solving linear system. We will use singular value 
decomposition (SVD) in this paper. 



3 SYNTHETIC LIBRARY 

We have computed a grid of synthetic spectra covering the 
wavelength range 354-716 nm with wavelength steps equiv- 
alent to 0.6 km s _1 . The calculations are based on Kurucz 
ODFNEW model atmospheres (Castelli & Kurucz 2004) and 
the synthesis code ASSeT (Koesterke 2009; Koesterke, Al- 
lende Prieto & Lambert 2009), operated in ID mode. The 
reference solar abundances for the synthesis are from As- 
plund, Grevesse & Sauval (2005), and while temperature and 
density are taken from the model atmospheres, the electron 
density is recalculated for consistency with the equation of 



'- 



1 r 



V : 

oi 3 L 

Cn : 

O : 

4 1 



5r ~ -E 

30000 20000 14000 9500 6400 4200 3000 
Teff (K) 

Figure 1. Coverage in surface temperature and gravity of the 
model spectral grids presented in this work. For each combination 
of these two parameters there are 7 models covering the metallic- 
ity range between —2.5 and 4-0.5 dex with steps of 0.5 dex. 



Table 1. Range of parameters for each subgrid. 



Grid 


Teff 


log g 


[Fe/H] 


MILES stars 


# 


(K) 


(cm s — 2 ) 


(dex) 


included 


1 


3500:4500 


0:5 


-2.5:+0.5 


190 


2 


4750:6000 


0:5 


-2.5:+0.5 


350 


3 


6250:8000 


1:5 


-2.5:+0.5 


124 


4 


8500:11500 


2:5 


-2.5:+0.5 


48 


5 


12000:26000 


3:5 


-2.5:+0.5 


22 



state used, which includes the first 92 elements in the pe- 
riod table and 338 molecules (Tsuji 1964, 1973, with some 
updates). Partition functions are adopted from Irwin (1981). 

Bound-free absorption from H, H~, Hel, Hell, and the 
first two ionisation stages of C, N, O, Na, Mg, Al, Si, Ca 
(from the Opacity Project; see, e.g., Cunto et al. 1993) and 
Fe (from the Iron Project; Bautista 1997; Nahar 1995) is in- 
cluded. Line absorption is included in detail from the atomic 
and molecular (H 2 , CH, C 2 , CN, CO, NH, OH, MgH, SiH, 
and SiO) files compiled by Kurucfl The radiative transfer 
calculations account for Rayleigh (H) and electron scatter- 
ing. We considered stars in the ranges 3500 ^ T c s ^ 26, 000 
K, loggr < 5, and -2.5 <[Fe/H]^ +0.5, divided into 5 
subgrids, as shown in Table 1 and Figure 1. 

The spectra were smoothed to a resolution of 0.23 nm, 
and resampled in steps of 0.09 nm, matching the resolution 
and wavelength scale of the MILES library. 



4 CALIBRATION RESULTS 

We apply the fitting procedure described in ij2] to model 
the ratio of the observed to the synthetic spectra with low 
order polynomials. The procedure is independently applied 
to each of the subgrids in Table 1. We tested with various 
combinations of first, second, and third order polynomials 

3 kurucz.harvard.edu 
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Figure 2. Comparison between the atmospheric parameters for the MILES stars in the range of grid # 2 (see Table [3). The input 
parameters are those provided with the library (Cenarro et al. 2007), and the output parameters those rederived in this work. In the 
upper panels, the solid and dashed lines are straight lines with slopes unity and determined from least-squares, respectively. The bottom 
panels show histograms of the residuals, and Gaussians with parameters derived again from least-squares fitting. 



for each parameter, finding that second-order or higher de- 
pendencies for some parameters, even if improved agreement 
with the observations for the stars under consideration, led 
to unphysical shapes in poorly constrained regions of the 
parameter space. 

After some experimentation, the adopted polynomials 
were second order in temperature (n = 2), zero-th order in 
surface gravity (m = 0), and up to first order in metallicity 
(h = 1 for grids #1, 2 and 3, but h = for grids #4 and 
5). The polynomial model for the ratio of observed and syn- 
thetic spectra was then used to correct the synthetic grid, 
and new synthetic spectra for the parameters of the MILES 
stars were derived by quadratic Bezier interpolation. As ex- 
pected the corrections tightened, made more symmetric, and 
centred closer to zero the distributions of residuals. 

Particularly large residuals remain after the correction 
for some stars with cool temperatures and low gravities, sug- 
gesting that these stars are somewhat singular, or that their 
assigned parameters are wrong. To improve consistency, we 
redetermined the atmospheric parameters for all stars by 
using the most recent version of the optimisation code dis- 
cussed by Allende Prieto et al. (2004, 2006, 2008, 2009), and 
exclude a few stars which could not be fit reasonably well. 
Fig. 2 illustrates the comparison between the old and the 
rederived parameters for grid # 2. Robust estimates of the 
mean and the standard deviation (half of the width of the 
distribution after discarding 15.85 % of the sample on each 
end) between the MILES parameters and the redetermina- 
tions are given in Table 2, which also includes the numbers 
of surviving library stars within each of the subgrids. 



Table 2. Standard deviation between the original and updated par. 



Grid 


T e ff 


logg 


[Fe/H] 


MILES stars 




(K) 


(cm s 2 ) 


(dex) 


surviving 


1 


94 


0.48 


0.26 


166 


2 


98 


0.31 


0.19 


337 


3 


128 


0.19 


0.17 


118 


4 


616 


0.30 


0.51 


45 


5 


783 


0.19 


0.39 


19 



The parameters provided with MILES (Cenarro et al. 
2007) have been compiled from the literature, and subse- 
quently homogenised by identifying and removing systemat- 
ics across data sets (Cenarro et al. 2001). The metallicities 
in this compilation are mainly from high-resolution stud- 
ies, and therefore reflect measurements from iron lines. Our 
model spectra have simply solar scaled metal abundances 
and, interestingly, if the corrections introduced are indepen- 
dent of [Fe/H], i.e. h — 0, our rederived metallicities are 
systematically higher than those in the MILES catalogue for 
metal-poor stars. With h > 0, such a trend disappears, as 
the corrections partially account for the relative strength- 
ening of the features produced by a— elements relative to 
iron. 

The polynomial fitting using the algorithm described 
in Section [5] was then repeated using the original grid of 
synthetic spectra and the rederived atmospheric parameters. 
New model fluxes for the MILES stars were recalculated by 
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Figure 3. Residuals for all the stars in each subgrid, and all wave- 
lengths, using the original synthetic libraries (solid black line), the 
corrections using n = m = h = 1 (solid dashed line), those for 
n = 2, m = 0, h = I (red), and those for n = 2, m = 0, h = 
(blue). 
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Figure 4. The lines show the 25% and 75 % percentiles of the dec- 
imal logarithm of the absolute value of the residuals between the 
model and the observed spectra for the original theoretical library 
(black) and the 2nd-order corrected (lst-order for the warmest 
subgrid, shown in the bottom panel) library (red). 



quadratic Bezier interpolation, and the residuals between 
these and the observed fluxes are displayed in Fig. 3. 

Fig. 4 shows the 25% and the 75% percentiles of the 
absolute value of the residuals for each of the subgrids as a 
function of wavelength. The black lines are the original re- 
sults for the theoretical library, and the red lines correspond 
to the same percentiles for the corrected library. Not only 
the residuals are decreased for essentially all the parameter 
space, but they are also smoother with wavelength for cool 
spectral types. 

Fig. [5] illustrates the same results from a different per- 
spective. The absolute value of the residuals are now shown 
for individual stars for each subgrid. Each star has two data 
points in the plot, at the same location, differing only on 
the size and type of symbol: open circles correspond to the 
original synthetic grid, and filled circles to the second-order 
corrected grid. The corrected library is in generally closer to 
the observations, but not in all cases, as permitted by the 
least-squares approach. 

To ensure that we are not overfitting the data, the poly- 
nomial modelling was repeated using only 2/3 of the stars, 
examining the changes in the residuals between the corrected 
grid and the observations for the remainder of the stars. This 
test showed also a significant reduction in the residuals, in 



line with the results obtained when fitting the complete sam- 
ple. 



5 EXTERNAL CHECK 

A sanity check of the empirical spectral energy distributions 
can be done by using stars with reliable spectrophotometric 
calibrations. We have chosen three key stars: the Sun, BD 
+ 17 4718, and Vega. Our choice of atmospheric parame- 
ters for these stars, compiled from the literature, are shown 
in Table 3. The data are from the calibration archive for 
HST, calspec (Bohlin 2008, 2010, and references therein); 
they come from HST observations for all the stars but the 
Sun, which are from the compilation by Colina, Bohlin, & 
Castelli (1996). We have dereddened the spectrum of BD 
+ 17 4718 according to the prescription by Fitzpatrick & 
Masa (Fitzpatrick 1999), using A V /E(B - V) = 3.1, consis- 
tent with the corrections applied to MILES. 

We warn the readers that parameters adopted for BD 
+ 17 4708 and Vega rely on the HST spectrophotometry 
for these stars and model fluxes based on Kurucz atmo- 
spheres, and therefore the excellent agreement found for our 
purely theoretical fluxes is somewhat facilitated. The test is 
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Figure 5. Absolute value of the residuals for individual stars in each of the grids. The open circles correspond to the original grid of 
synthetic spectra, and the filled circles to the grid corrected with first-order polynomials. The residuals are projected onto the T e fj-log<; 
plane for all five grids, and also for the T e g-[Fe/H] plane for grid #2 in the bottom-left panel. 



Table 3. Testing stars. 



Star 


T c ff 


log9 


[Fe/H] 


E(B-V) 


References 




(K) 


(cm s 2 ) 


(dex) 


(mag) 




Sun 


5777 


4.44 


0.0 


0.00 


Stix (2004) 


BD +174708 


6141 


3.87 


-1.7 


0.01 


Ramirez et al. (2006) 


Vega 


9620 


3.98 


-0.7 


0.00 


Garci'a-Gil et al. (2005) 



only completely fair for the Sun, as this is the only star in 
this group with parameters that are truly independent from 
model atmospheres (and have negligible errors), but the oth- 
ers are included as they may help to confirm the conclusions 
found from the examination of the solar case. 



Fig. 6 compares the predicted spectral energy distribu- 
tions, obtained by linear interpolation in a) our purely the- 
oretical grid of model spectra (red line in left-hand panels), 
and b) our empirically calibrated grid (red line in right-hand 
panels). The interpolated spectra have also been smoothed 
to approximate the lower resolution of the HST spectropho- 
tometry, with a FWHM resolving power A/<5A ~ 1000. Over- 
all, the calibrated grid performs slightly better than the 



purely theoretical grid in some regions, but the opposite is 
true in others. 



This result suggests that there may be some systematic 
differences between the flux scale adopted for HST calibra- 
tion and that of the MILES library, and the corrections de- 
rived from MILES may not be widely applicable. To keep a 
perspective, we note that systematic and random errors in 
the MILES fluxes were estimated to be about 2 and 3 %, 
respectively, from the comparison with photometry (B — V) 
from the Lausanne data base (Mermilliod et al. 1997). A re- 
duction in effective temperature for the solar model by 150 
K - which is a typical error in this quantity to be expected 
for most stars (not the Sun)- would be enough to bridge the 
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5% gap in the red between the solar fluxes and those from 
the empirically corrected library. 

It should be emphasised that an empirical calibration, 
such as the procedure described in this paper, will equally 
erase discrepancies between models and observations due to 
model deficiencies and due to systematic errors in the obser- 
vations, and the latter are highly dependent on the source of 
the data. In the upper-right panel of Fig. [6l we also show in 
blue the spectra of three G2V stars included in the MILES 
library: HD 10307, HD 84737, and HD 13043. As expected, 
their fluxes are consistent with the spectrum interpolated in 
the corrected librarjQ. An interpolated spectrum for the so- 
lar parameters from neighbouring stars using the tool avail- 
able at the MILES web sit^fl leads to the same conclusion. 

It is interesting to highlight the poor matching between 
the original library and the observed strength of the Ca II 
H and K lines in BD +17 4708. This mismatch, which can 
be directly ascribed to the use of solar abundance ratios 
in the models, and hence the neglect of the enhancement 
in a elements typically observed in metal-poor halo stars, 
is reduced somewhat after the polynomial corrections have 
been introduced. 



6 SUMMARY AND CONCLUSIONS 

We present a general strategy to calibrate a set of model stel- 
lar fluxes using observations for a network a flux standards. 
At each frequency, the ratio of model and observed fluxes 
is fit by least-squares to a polynomial that depends on the 
atmospheric parameters. We test this scheme on the MILES 
stellar library, with satisfactory results: second-order cor- 
rections on effective temperature, zero-th order corrections 
on surface gravity, and zero-th or first-order corrections on 
metallicity perform well and improve significantly the agree- 
ment between the model grid and the observations, resulting 
in tighter and more symmetric distributions of residuals. 

Seeking an external assessment of the performance of 
the corrections, we take a close look at the spectra of three 
well-known flux standards (see Bohlin 2007 and references 
therein). The result of this test is a much modest than ex- 
pected improvement, showing that the purely theoretical 
fluxes perform better than those corrected for some wave- 
lengths and stars. This suggests that the HST fluxes and 
those in the MILES library are not fully compatible, but 
further investigation is warranted. 

The HST/STIS Next Generation Spectral Library 
(Gregg et al. 2005) has been recently released as part of 
the high-level products in the Multi-mission Archive at the 
Space Telescope (MAST). This library is based on obser- 
vations made in cycles 10 through 13 and includes nearly 
400 stars. This number is smaller than those in the MILES 
library and the resolution of the spectra is also lower, but 
the spectral coverage is wider. An application of the method 
presented here to this library, or if compatible with MILES, 
to the combination of the two libraries would be useful and 
is planned for the future. 

4 A fourth G2V star HD 76151 in the library, however, shows 
fluxes that more closer to the reference solar observations. 

5 http:/ /www.iac.es/proyecto/miles/pages/webtools/star-by- 
parameters . php 
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Figure 6. Observed (CALSPEC) spectrophotometry for three 
flux standards (Sun, BD +174707, and Vega; in black) compared 
with spectra interpolated from our grids (red lines): left-hand side 
graphs for the purely theoretical grid, and right-hand side for 
the MILES-calibrated grid. In the top right-hand panel, several 
spectra of solar-like stars from the MILES library are shown in 
blue. 
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